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ABSTRACT 

Observations performed in the last few years indicate that most massive early-type 
galaxies (ETGs) observed at redshift z > 1 exhibit sizes smaller by a factor of a few 
than local ETGs of analogous stellar mass. We present numerical simulations of the 
effect of baryonic mass loss on the structure of a spheroidal stellar system, embedded 
in a dark matter halo. This process, invoked as a possible explanation of the observed 
size increase of ETGs since z ~ 2, could be caused either by QSO/starburst driven 
galactic winds, promptly ejecting from Early Type Galaxies (ETGs) the residual gas 
and halting star formation (galactic winds) , or by stellar mass returned to the ISM in 
the final stages of stellar evolution. Indeed, we find that a conceivable loss of ~ 50% 
of the baryonic mass can produce a significant size increase. However, the puffing up 
due to galactic winds occurs when the stellar populations are much younger than the 
estimated ages £ 0.5 Gyr of compact high-z ETGs. Therefore, while it may have 

had a role in deciding the final structure of ETGs, it cannot explain the evolution 
observed so far of their size-mass relation; its signature should be searched for in 
much younger systems. Conversely, the mass loss due to stellar evolution could cause a 
relatively modest expansion of passively evolving stellar systems later on, contributing 
to, without dominating, the observed evolution of their mass-size relationship. 

Key words: galaxies: formation - galaxies: evolution - galaxies: elliptical - quasars: 
general - method: numerical 



1 INTRODUCTION 

During the last years it has been established that most mas- 
sive early-type galaxies (ETGs) observed at redshift z > 1 
exhibit sizes smaller by a factor of a few than local ETGs 
of analogous stellar mass (e.g. Daddi et al. 2004; Trujillo et 
al. 2006, 2007, 2011; Longhetti et al. 2007; Toft et al. 2007; 
Zirm et al. 2007; van der Wei et al. 2008; van Dokkum et al. 
2008; Cimatti et al. 2008; Buitrago et al. 2008; Damjanov et 
al. 2009; Ryan et al. 2010). 

The possibility that the size evolution is, at least in 
part, apparent and due to some subtle systematic effect has 
not been completely ruled out. Discussed caveats include a 
centrally concentrated source (such as an AGN or central 
starburst), age gradients, un-detected low surface bright - 
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ness external regions at high z, or a top-heavy IMF affecting 
mass estimates (see Daddi et al. 2005; Van Dokkum et al. 
2008; La Barbera et al. 2009; Mancini et al. 2010; Hopkins 
et al. 2009, 2010). Moreover, very recent observational re- 
sults claim some co-existence of compact and normal size 
ETGs, both at high (Mancini et al. 2010; Onodera et al. 
2010; Saracco et al. 2010) and low redshift (Valentinuzzi et 
al. 2010). 

However, nowadays most authors agree that the obser- 
vational results are dominated by a real size evolution. The 
proposed interpretations are related either to the effects of 
mergers or to the loss of a substantial fraction of mass from 
the galaxy. 

Khochfar & Silk (2006) presented a semi-analytic model 
where the size evolution is substantial only for galaxies more 
massive than 10 11 M©, and results from massive galaxies at 
high redshifts forming in gas-rich dissipative mergers, whilst 
galaxies of the same mass at low redshifts form from gas- 
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poor mergers. However, current observations indicate size- 
able evolution also in much smaller systems (e.g. Ryan et 
al 2010). Similar ideas have been explored by Hopkins et 
al. (2009) , in a more phenomenological model incorporating 
results of a large suite of numerical simulations of mergers. 

Actually, the size increasing effects of major (wherein 
the two merging galaxies have comparable masses) dry (i.e. 
without a significant collisional gas component) mergers has 
been often discussed. By converse, in wet mergers, the pres- 
ence of a dissipative gas component limits the gain in size. 
However, even the former process faces a few problems. The 
most basic is that it would move galaxies too slowly toward 
the local size-mass relationship (Cimatti et al 2008; Bezan- 
son et al. 2009; Damjanov et al. 2009; Saracco et al 2009). 
This is because, according to simple virial theorem argu- 
ment, confirmed by a number of numerical simulations, in 
major dry mergers the size increases almost linearly with the 
mass, and possibly somewhat less (e.g. Ciotti & van Albada 
2001; Nipoti et al. 2003; Boylan-Kolchin et al. 2006; Naab 
et al. 2009). This dependence is too close to the observed 
local mass-size relationship r oc M 0,56 (Shen et al. 2003) 
to explain the evolution in a reasonable number of merging 
events. 

The most promising merging mechanism to explain the 
size increase seems to be a series of late minor dry merging 
events. These would add stars in the outer parts of passive 
high-z galaxies, in such a way to produce a size increase that 
can scale as steep as M 2 (e.g. Naab et al. 2009; Oser et al. 
2010). Actually, Hopkins et al. (2010), considering the vari- 
ous proposed channels for observed size evolution, concluded 
that minor dry mergers are the best candidate to dominate, 
though other channels should have a non negligible role. 

In this paper, we will test the possible contribution of 
the puffing-up process. This envisages that the expansion in 
size is driven by the expulsion of a substantial fraction of 
the gas out of the galaxy either by AGN and/or supernova 
driven galactic winds (Fan et al. 2008, 2010), or by the expul- 
sion of gas associated to stellar evolution (Damjanov et al. 
2009). In the former case the expulsion timescale would be 
short, likely not much longer than the dynamical timescale, 
at least when driven by the AGN, whilst in the latter an im- 
portant mass loss could last even ~ 0.5 — 1 Gyr, depending 
on the IMF and on stellar evolution details. 

While the works by Fan et al. (2008, 2010) used as a 
reference the specific semi-analytic model by Granato et al. 
(2004) for the co-evolution of ETGs and SMBH (henceforth 
G04) , this kind of puffing-up due to baryonic mass loss from 
the galaxy has a broader applicability. Indeed, virtually all 
modern models of galaxy formation give a prominent role 
to AGN and/or SNae driven galactic wind, ejecting from 
the galactic region a substantial fraction of gas (for a recent 
review, see Benson 2010). Thus, it is very likely that this 
puffing up played a role in deciding the final sizes of ETGs, 
at some point over their history. However, it is still an open 
question whether this role has been major, and, in particu- 
lar, whether it can explain the available observations. The 
aim of the present work is to provide a step to clarify it. 

The idea that a puffing-up mechanism should at some 
point have had occurred rests on many hints, nowadays com- 
ing from a complex interplay between observations and the- 
ory, suggesting that ETGs are the result of an intense phase 
of high-z star formation activity, possibly traced by the sub- 



mm selected galaxy population. It is unlikely that these ex- 
treme "starburst", induced by fast collapse and/or by rapid 
merging of gas rich systems, have been suddenly terminated 
by simple gas consumption. It is instead usually envisaged 
that some strong feedback was capable to eject from the 
galaxy, on a relatively short timescale, a substantial fraction 
of its baryonic mass, in the form of gas not yet converted 
into stars (a process usually referred to as "galactic wind 
or superwind"; see e.g. Benson et al. 2003; Pipino, Silk & 
Matteucci 2009; for observational indications of these high- 
z galactic outflows see Lfpari et al. 2009 and Prochaska & 
Hennawy 2009, and references therein). In the past years it 
has been pointed out that the most likely candidate to power 
such a process, at least in massive systems, is QSO activity 
(e.g. Silk & Rees 1998; Fabian 1999; Granato et al. 2001, 
2004; Benson et al 2003; Cattaneo et al. 2006; Monaco et al. 
2007; Sijacki et al. 2007; Somerville et al. 2008; Johansson 
et al. 2009; Ciotti, Ostriker & Proga 2009). 

Additionally, after the termination of this huge star- 
forming phase, it is likely that the galaxy lost another sig- 
nificant fraction of its baryonic mass, due to stellar evolution 
(supernova explosions and stellar winds). 

Independently of the still uncertain details of the forma- 
tion mechanism of ETGs, or more generally of the spheroidal 
component of galaxies, it seems timely to investigate with 
aimed numerical simulations the effects of the expulsion of a 
significant fraction of baryonic mass from a spheroidal sys- 
tem embedded in a dark matter (DM) halo. 

The plan of the paper is the following. In Section [2] we 
recall previous results on a problem sharing similarities with 
that considered here: the dynamical evolution of young star 
clusters after dispersion of the parent gas cloud. These re- 
sults inspired the proposal that the size evolution of ETGs 
could be due to the puffing up mechanism, and were adopted 
as a rough approximation for quantitative considerations. 
Section [3] we describe the simulation technique and the ini- 
tial conditions, the results are presented in Section 2] and 
discussed, in the context of observed size evolution of ETGs, 
in the final Section [5] We use the concordance cosmology 
(Komatsu et al. 2009), i.e., a flat universe with matter den- 
sity parameter Qm = 0.3 and Hubble constant Ho — 70 km 
s _1 Mpc" 1 . 



2 THE STAR CLUSTER APPROXIMATION 

A process somewhat similar to the puffing up of ETGs by 
mass loss has been addressed several times, both analytically 
and numerically, in the context of the dynamical evolution 
of star clusters. Again, the mass loss is due to (i) the in- 
direct and combined effects of stellar winds and supernova 
explosions, soon driving out the significant fraction of gas 
not used in star formation, or (ii) directly to stellar evolu- 
tion. As for star clusters, the latter contribution includes the 
gas ejected from the stars by winds and explosions, but also 
the compact renmants that may get at birth a kick velocity 
sufficient to be ejected from the system. The most obvious 
differences, with respect to the puffing-up of ETGs, are the 
absence of an embedding dark matter halo and, to a lesser 
extent, the importance of two-body collisions. 

Biermann & Shapiro (1979) and Hills (1980), under sim- 
plifying assumptions to allow an analytical treatment, found 
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that the size evolution depends on the ejection timescale. If 
this is short compared to the dynamical time (hereafter fast 
ejection), conservation of specific kinetic energy yields for 
the expansion factor 

Hi 



(1) 



2e- 1 

in terms of the ratio between the final and initial mass 
e = Mi/Mi; note that for e ^ 0.5 the system becomes un- 
bound and dissociates. On the other hand, if the ejection 
timescales is much longer than the dynamical one, conser- 
vation of adiabatic invariants yields the size evolution 

Rt_ _ 1 

Ri ~ e ; 



(2) 



Thus, fast expulsion is more effective in increasing the size, 
while, if the expulsion is slow enough (adiabatic), the system 
remains bound independently of e. 

These relationships have been checked and substantially 
confirmed by several numerical simulations of star clusters 
dynamics, starting from the pioneristic work by Tutukov 
(1978). However, due to a few effects not accounted for in 
analytical works, a portion of the system remains bound 
even if e is somewhat smaller than 0.5 and the mass loss is 
fast (see Baumgardt & Kroupa 2007 and references therein). 
Numerical experiments show also that, after fast mass loss 
ends, the system rapidly reaches a maximum transient ex- 
pansion within 10 — 15 dynamical times, while a new equi- 
librium is attained over 30 — 40 dynamical times (e.g. Geyer 
& Bukert, 2001; Goodwin & Bastian, 2006; Baumgardt & 
Kroupa, 2007). 

Fan et al. (2008, 2010) tentatively tested against avail- 
able data a QSO driven puffing up scenario for ETGs, resting 
on the G04 model of joint evolution of SMBH and spheroids, 
and adopting the above relationships coming from star clus- 
ter dynamics. Their findings are to some extent encourag- 
ing, but, as remarked by Mancini et al. (2010), a closer 
analysis reveals that the expansion time-scale seems far too 
short to explain the relatively old ( ^ 1 Gyr) stellar ages 
claimed for high-z compact galaxies. A similar, albeit less 
dramatic, problem has been pointed out by Damjanov et al. 
(2009) when trying to explain the observed size evolution by 
means of mass loss due to stellar evolution. In this case the 
timescale of mass loss are dictated by stellar evolution, and 
the expected expansion is adiabatic. 

However, assuming the above recipes for ETGs is only 
a zeroth order approximation, since in ETGs the DM halo 
is expected to affect the efficiency and timescale of the size 
evolution, and to prevent galaxy disruption even when a 
major fraction of baryonic mass is lost. In the following, our 
purpose is to investigate these effects via simple but aimed 
numerical simulations. 



3 NUMERICAL METHOD AND SETUP 

The purpose of the simulations is to investigate the evolu- 
tion of collision-less particles (stars and DM) under a change 
of gravitational potential due to a loss of baryonic mass of 
the system. The escaping mass can be either the gas which 
has not been converted into stars during the star forming 
phase of the spheroid, or the mass lost from stars in form of 
stellar winds and SNae explosions. In any case, we assume as 
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Figure 1. Half mass radius increase R{/Ri due to baryonic mass 
loss as a function of the parameter e = Mgji n /Mg ini , for dif- 
ferent ejection times At (thick lines are with DM included). The 
thin lines illustrate the size increase according to the analytic ap- 
proximation, meant for a system not embedded in a DM halo, of 
Eq. (1) (solid, fast ejection) and (2) (dashed slow ejection). The 
points show examples of runs without DM halo included (aster- 
isks for fast ejection At = and circles for slow ejection At = 80 
Myr) 



given, and due to "external" causes (such as SNae and AGN 
feedbacks, or stellar evolution), the temporal dependence of 
this mass loss fEq, ll2[) , which we put by hand, and we simu- 
late the ensuing evolution of collision-less mass distributions. 
Therefore we don't have to treat the gas dynamics. This is 
the same approach followed in most simulations of puffing- 
up of star clusters (e.g. Boily & Kroupa, 2003; Goodwin & 
Bastian, 2006; Baumgardt & Kroupa, 2007). 

We used the A^-body code GADGET-2 (Springel et al. 
2005) to perform simulations with 10 6 and 5 x 10 6 particles. 
None of the presented results show any noticeable difference 
in the two cases, which assures us that the mass resolution 
is sufficient for the purposes of the present study. Half of the 
particles are used to sample the baryonic and dark matter 
components respectively. 

The density distribution of DM particles is assumed to 
follow the standard NFW (Navarro, Frenk & White 1997) 
shape 



Pdm(t) = 



Mvir 



! g(c) 



47r/?J ir r(l+cf) 



(3) 



where M v j r ,DM is the halo virial mass in DM (the DM mass 
inside ii v j r ), f — r/R vlI , c is the concentration parameter 
and g{c) = [log(l + c) - c/(l + c)]" 1 . 

The virial radius, i? v ir, is by definition the radius within 
which the mean density is A V i T (z) (a quantity classically 
coming from the spherical top-hat collapse model) times the 
mean matter density of the universe p u (z): 
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Figure 2. Ratio R/Ri of the current to the initial half-mass radius as a function of time, for different values of the diet parameter e 
and of the ejection times At, as indicated in the panels. The thin solid line in the top right panel shows the evolution of a model not 
including the DM halo, for e = 0.6. Note that the latter is shown in the panel corresponding to e = 0.4 (for the cases including DM) since 
the vertical scale is more adequate. The double arrows show the duration of a 10 td yn time interval. For the adopted initial conditions, 
t dyn ~ 5 Myr. See text for explanations. 



Rvir 



47T A v ir(z)p u (z) 



1/3 



(4) 



The overdensity A v ir(z), for a flat cosmology, can be 
approximated by 



A vir (2:) 



(187r 2 +82x- 39a; 2 ) 



(•») 



where x — Q(z) — 1 and Q(z) is the ratio of the mean 
matter density to the critical density at redshift z (Bryan & 
Norman 1998). 

The corresponding mass distribution is written 



M D m(< r) = M v ir,DM3(c) 



log (1 + cf) 



(6) 



(1 + cf)_ 

To cope with the divergence of the NFW mass distribution, 
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we introduce an exponential cut at 3i? v ir, which very safely 
contains the region where the baryonic component is im- 
portant. We performed also a few test runs with the cut 
at 2i? v i r , verifying that none of the discussed quantities is 
affected by this choice. 

For the baryonic particles (stars and gas), we assume 
that they follow an Hernquist (1990) profile, which provides 
a reasonable description of stellar density in local spheroids: 



, . Mb a 



(7) 



2ir r (r + a) 3 ' 
The corresponding mass distribution is 

M B (<r) = M B Y ; (8) 

so that the half-mass radius is related to the scale radius a 
by R1/2 = (1 + \/2) a and, assuming a mass to light ratio 
independent of r, the effective radius is R e ~ 1.81a. 

In the following, unless otherwise specified, by dynam- 
ical time idyn we mean the initial (i.e. before any mass loss 
and expansion) dynamical time, computed at Ri/2- 



idyn 



2G (M b /2 + Mdm(< Ri/2)) 



1/2 



(9) 



For our standard initial conditions (see below) the contri- 
bution of DM to the mass inside R1/2 amounts to ~ 20%. 
Thus idyn can be estimated (within 10%) neglecting it: 



t-dy 



2.3 



f-^Y 

Vikpcy 



Mb 



10 n M< 



Myr 



(10) 



Given the density runs, we obtain the ID velocity dis- 
persion by integrating the Jeans equation under the assump- 
tion of isotropic conditions: 



Px(r) 



' ,GM TOT {<r') , 
dr ^ px{r ) 



(11) 



where X stands for B or DM, and Mtot(< r) = Mdm(< 
r)+Ms(< r). By evolving the particle system for several dy- 
namical times, we get confident that it is actually in (quasi- 
static statistical equilibrium. 

Starting from this initial setup, we introduce a mass 
loss, intended to emulate the various possible effects de- 
scribed above, by removing exponentially over an ejection 
time At a fraction 1 — e of the baryonic mass : 



M B (t) = M B ( t =o) exp 



(£') 



(12) 



For instance, this simple functional form provides an accept- 
able description of the gas removal due to QSO feedback in 
the G04 semi-analytic model, with an ejection timescale of 
the order of 20-30 Myr for a wide range of the model pa- 
rameters. 

The mass loss is practically attained by decreasing cor- 
respondingly in time the mass of the baryonic particles sam- 
pling the density field. After the end of the mass loss period, 
we let the system to evolve till it reaches, if any, a new equi- 
librium configuration. 

The reference value for the initial (i.e. before any mass 
loss) ratio of virial mass (total mass within the virial radius) 
to baryonic mass is M v i r /M B ( t =o) — 25. In a recent analy- 
sis consistent with previous works, Monster et al. 2010 find 



that this ratio should be, in the local universe, about 50 for 
DM haloes of ~ 5 x 10 12 M©. However, to get a significant 
puffing-up, the system should lose something of the order of 
50% of its baryonic mass. Thus, we set as initial condition 
a ratio twice smaller than that. 

We set M v i r = 10 13 Mq in all simulations. Nevertheless, 
our results apply to different values of M v i r , provided that 
the ratios of scale radii and masses in the two components 
(DM and baryons) are not changed, and that the time is 
measured in units of dynamical time idyn oc p^ 1 ^ 2 . 

We adopt a concentration parameter c = 4, a typi- 
cal value at galactic halo formation (see Zhao et al. 2003; 
Klypin, Trujillo- Gomez, & Primack 2010), and i? vir ~ 170 
kpc, which is the value given by Eq. [4] and Eq. [5] for a 
M vir = 1O 13 M halo virialized at z = 3. 

We set a = 1.5 kpc (Re ~ 2.7 kpc). This seems a value 
adequate to study the evolution of the system in the plane 
effective radius R e vs. stellar mass M*. Indeed, assuming 
that about half of the initial baryonic mass is in form of 
stars, the system would lie initially a factor ~ 2.5 below 
the local mass-size relationship for ETGs. The initial (i.e. 
before mass loss and expansion, Eq. |§J dynamical time is 
idyn ~ 5 Myr. Note that a smaller initial size would shorten 
the dynamical time, exacerbating the problems pointed out 
in Section [5] 

In summary, the parameters affecting the results of our 
simulations are the ratio of mass between the total and 
baryonic components M v i r /M B ( t=0 ); the corresponding ra- 
tio of scale- lengths R vlI /a; the fraction of baryonic mass 
lost (1 — e), and the time At over which the loss occurs. We 
performed simulations covering broad ranges of the latter 
two quantities, while in most runs we kept the former two 
at the fiducial values discussed above. We checked however 
that none of our qualitative conclusion is affected by factor 
~ 2 variations of them, and likely even by larger ones (see 
discussion at the end of Section [4}. 



4 RESULTS 

As a preliminary sanity check, we ran simulations without 
DM with ejection times At = Myr and 80 Myr, i.e., 
much shorter or much longer than td yn , respectively. The 
ratio of the final to the initial half-mass size Ri/Ri gen- 
erally agrees well with the expectations given by Eq. (1) 
and (2) (see Fig. [TJ. However, when the mass loss is fast 
and approaches 50%, the analytical formula Eq. (1) increas- 
ingly over-predicts the numerical result, underlying the well 
known fact that the divergence does not occur at e = 0.5 
in numerical simulations. This is in keeping with previous 
findings (e.g. Geyer & Burkert 2001). 

We now turn to cases with DM included. Fig. [1] illus- 
trates the size expansion as a function of the fraction of 
remaining baryonic mass e = 0.2, 0.4, 0.6, and 0.8, for dif- 
ferent ejection times At = 0, 2, 5, 20, and 80 Myr. The ex- 
pansion increases with decreasing e and At, but it is milder 
with respect to the corresponding purely baryonic case. In 
particular, the system is no longer disrupted even when the 
ejection is impulsive (At = 0) and e is as low as =0.2; this is 
expected since DM constitutes the dominant source of gravi- 
tational potential at large radii. We have verified, with a few 
sample runs, that the case with At = 80 Myr is represen- 
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Figure 3. The time needed, after the end of mass loss, to re- 
cover an equilibrium configuration as a function of the parameter 
e = M B j in /Mg i„i, for different ejection times At. This time 
is defined by the epoch after which the size never changes more 
than 10% with respect to the final value. 



tative of the "slow expulsion" regime At 3> tdyn- In other 
words, the expansion factor do not decrease any more for 
larger values of At. Note also that, as already mentioned, 
At — 20 is likely the case that best approximates, among 
those shown, the QSO driven gas expulsion predicted by the 
G04 model (see Fig. 0), and considered by Fan et al. (2008; 
2010). The corresponding expansion is significantly smaller 
than that achieved for instantaneous ejection. 

Fig. [5] shows in detail the time evolution of the system 
size during and after the ejection. For sake of comparison we 
include also a case without DM. For At 3> tdyn, the size in- 
creases mostly during the mass loss, and stabilizes soon after 
At, so that the system is actually in quasi-equilibrium dur- 
ing the ejection process; contrariwise, for impulsive ejection 
with At <C tdyn the initial equilibrium is totally broken, and 
the system expands significantly more. The size undergoes 
damped oscillations before stabilizing, more important for 
smaller e or At. Fig. shows the time needed, after the end 
of mass loss, to recover a substantially stable configuration. 
This turns out to be shorter than the corresponding time 
without the DM component, as expected on an intuitive ba- 
sis, due to the stabilizing effect of the latter (for instance, 
compare the thin solid curve in the top-right panel of Fig. 
[2 with the thick solid curve in the bottom-left panel). 

Fig. [4] shows the evolution of the average stellar veloc- 
ity dispersion. During the size expansion the system cools 
down and the stellar random motions slow, reducing the av- 
erage dispersion. The net effect is found to be much evident 
for strong ejection (small e) but almost independent of the 
ejection time At. 

In Fig [5] we illustrate the effects of baryon expulsion on 
the density profiles of both components. For e 0.6 and 



fast expulsion, the baryonic component, after a violently 
disturbed phase, eventually recovers a density distribution 
reasonably well described by the Hernquist formula, albeit 
with a larger scalelength (for instance a = 2.34 with reduced 
X 2 = 2.6 for e = 0.6, shown in the figure). By converse, at 
lower e, the Hernquist fit becomes increasingly unaccept- 
able for the final equilibrium profile (e.g. a — 5 and reduced 
X 2 = 41 for e = 0.2). This is not surprising, since the fi- 
nal bound state is increasingly dictated by the presence of 
the embedding DM halo. Actually for e ,$ 0.5 the system 
would dissolve if not for the halo. In any case slower expul- 
sion with same e yields lower deviations from the Hernquist 
functional form. In particular, for e ^ 0.6 and At = 80 
Myr the Hernquist formula provides a good description of 
the density distribution even during the mass loss-expansion 
phase. 

For DM, on large scales the profile is unaffected, while in 
the inner region the baryon expansion drags an expansion of 
the DM particles. As a result, the DM profile in the galactic 
region is always flattened to some level with respect to the 
original NFW shape. 

Fig [6] shows the sensitivity of our results to the param- 
eters of the initial baryonic configuration. As expected, the 
effects are in the sense that, whenever the DM contributes 
more (less) to the mass inside the region occupied by the 
baryonic system, the latter expands less (more). This may 
occur by increasing (decreasing) the ratios M v i r /M B ( t = ) or 
a/Rvir. Note the trade-off between variations in initial size 
and its amplification due to mass loss. As a result, the final 
size is relatively insensitive to the initial one, particularly for 
fast expulsion. For instance, when the impulsive mass loss is 
60% (e = 0.4), a factor 4 change of the adopted initial R e , 
yields only a factor ~ 1.5 change in the final R e (blue points 
and arrows in Fig [6} ■ This could have a role in explaining 
the relatively low scatter of the observed local mass-size re- 
lationship. Also, an increase (decrease) of M B (t=o), keeping 
fixed all other parameters, and in particular M V i ri yields a 
larger (smaller) expansion. This tends to keep the expansion 
close to that required to maintain the final state of the sys- 
tem in the region of local mass-size relation (red points and 
arrows in Fig|SJ). 



5 DISCUSSION AND CONCLUSIONS 

An inspection to the previous figures, in particular to Fig. 
[2] (and to a lesser extent Fig. [3]), reveals the main prob- 
lem to explain the observed size evolution of ETGs with the 
puffing-up scenario. On one hand, our simulations confirm 
that, even in presence of a DM component, a factor ~ 2 in- 
crease in size can be expected in any galaxy formation model 
in which the spheroid loses ~ 50% of its baryonic mass, in 
particular when this happens on a timescale of the order of 
the dynamical timescale. Moreover, the expansion is larger 
for systems initially more compacts, resulting in an interest- 
ing self regulation of the final size (Fig. |6]). However, if this 
mass is constituted by the star forming gas, in scenarios in 
which a galactic wind suddenly sterilizes the galaxy (such 
as in the G04 model considered by Fan et al. 2008, 2010), 
the puffing-up occurs far too close to the last episode of star 
formation. Indeed, the galaxy is predicted to be smaller than 
the final size only for a very short time after expulsion, less 
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Figure 4. Left panel: Ratio o-f/ui of the final to initial mean 3D stellar velocity dispersion as a function of e, for different ejection times 
At. Right panel: same as a function of time, for different values of the diet parameter e (as labeled) and of the ejection times At. 




Figure 5. Two examples of the effects of baryon expulsion on the density profiles (shown as radial distribution of particles) of baryons 
(top panels) and DM (bottom panels, showing only the inner region). The left column refers to e = 0.6 and At = 0, while the right 
one is for e = 0.2 and At = 0. We plot three snapshots, i.e. the initial one, one wherein the baryons are far from equilibrium, and the 
final one (times indicated in the bottom panels). The blue solid lines are \ 2 fits to the distributions with Hernquist and NFW profiles 
respectively, which we show in the former case only for the initial and final states, and in the latter only for the inital one. 
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Figure 6. Initial position (points) of our runs in the stellar 
mass-size plane, and its evolution (arrows) due to gas ejec- 
tion amounting to 60% of the initial baryonic mass (i.e. for 
e = Mgji n /Mg i n i = 0.4). The various panels refer to different 
expulsion times At, as indicated. The starred point with black 
arrow represents the reference model, while the other four points- 
arrows show the behaviour of systems having an initial baryonic 
mass or radius twice smaller or greater (while the DM halo re- 
mains identical). For ease of reading, we have artificially displaced 
slightly in mass the points corresponding to variations in radius. 
The solid diagonal line is the mass-size relationship for local el- 
lipticals (Shen et al. 2003), with the 1 a dispersion depicted by 
dotted lines. The upper panel shows the lines above which the 
initial dynamical time is > 50 Myr, without DM (dot-dashed, 

Eg. 1101) or including it (dashed, Eq. [9jl (see Section 15) . 



than a few dynamical times, i.e. less than ~ 20 — 30 Myr 
for the adopted initial configuration. This is at least a fac- 
tor 20 less than the estimated ages of stellar populations in 
high-z compact galaxies ( Si 0.5— 1 Gyr; e.g. Longhetti et al 

2007; Damjanov et al. 2009). Even taking into account gen- 
erous uncertainties of these estimated ages, it seems clear 
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Figure 7. Example of evolution of the total baryonic mass (star 
forming gas+stars) in the G04 model for co-evolution of SMBH 
and spheroids (solid line, left axis), and the increase in size 
(dashed line, right axis) according to our corresponding simula- 
tion. The mass is measured in units of M max , which is the max- 
imum value of the total baryonic mass (stars+cold gas) attained 
during its formation. This is reached just before the galactic winds 
eject the cold gas and stop further star formation. The abrupt de- 
crease of the mass after ~ 0.3 Gyr marks the ejection of the cold 
gas, not yet converted into stars, by the AGN-driven wind, while 
the subsequent slow decrease is due to stellar mass returned to 
the ISM, under the assumption that the galaxy potential cannot 
retain it. See text for more details. 



that a substantial contribution of galactic winds to the size 
evolution observed so far can be safely ruled out. Neverthe- 
less, this process should have had a role in deciding the size 
of ETGs, if (QSO-driven) galactic winds caused their sud- 
den death, but its signature should be searched for in much 
younger systems. This poses huge observational challenges 
with present facilities. 

Since the expansion timescale is proportional to the 
dynamical time i<jyn oc M~ 5 R 15 , it could be suspected 
that the problem originates from our choice of initial con- 
ditions, such as an excessive initial baryonic mass or an in- 
sufficient initial radius. However to get a factor Si 10 in- 
crease (a minimal requirement) of the dynamical time, the 
initial size should be increased at least by a factor Si 5, or 
the mass decreased by a factor Si 100. Actually, even these 

large changes of the baryonic distribution would not suffice, 
since in both cases the fixed DM component would become 
important to set the dynamical time in the galactic region 
(Eq.[9]). Anyway, in both cases the initial state of the system 
would already lie well above the observed local size-mass re- 
lationship. This is illustrated by lines in the top panel of 
Fig. [6] delimitating the region of the plane R e — M* where 
^dyn Si 50 Myr. As can be appreciated there, the inclusion 
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of a DM contribution to estimate td yn makes the argument 
even stronger. Moreover, our standard initial conditions put 
the system about a factor 2.5 below the z=0 size-mass rela- 
tionship, if about half of the baryons are in stars (in other 
word if e ~ 0.5). A non negligible fraction of compact high-z 
ETGs are found up to a factor ~ 3 — 5 below this relation- 
ship (e.g. Zirm et al. 2007, Toft et al. 2007, Van Dokkum 
et al. 2008, Cimatti et al. 2008, McGrath et al. 2008, Cas- 
sata et al. 2010). For these, the initial dynamical time is 
even shorter, exacerbating the problem, as we anticipated in 
Section [3] On the other hand, given the anti-correlation be- 
tween initial size and expansion Rf/Ri achieved after mass 
loss, discussed at the end of Section [4] and shown in Fig. [(J 
also for these extremely compact objects it is conceivable to 
get final sizes close to the locally observed ones. 

We have also verified, with a few sample simulations, 
that this general conclusion on the shortness of the expan- 
sion timescale is not significantly affected by assumptions 
such as that the gas and stars share the same profile be- 
fore ejection of the former, or that this ejection occurs ho- 
mogenously in the system. The same holds true for different 
choices of the density profile of the DM and baryonic com- 
ponents (Eq. [3] and Eq. 0), provided they describe in a rea- 
sonable manner the density distributions of the respective 
mass distribution. 

Even in the case of expansion driven by stellar mass 
loss the problem of the excessive shortness of the expansion 
timescale is likely important, though less clear cut. Indeed, 
Damjanov et al. 2009 pointed out it, basing their reason- 
ing on the star cluster approximation reviewed in Section [2l 
and in the performing simplified estimates of mass loss due 
to stellar evolution. They found that the fraction of mass 
lost during the passive evolution of stellar populations can 
be as large as 30-50%, depending on the IMF, but the ma- 
jority of this loss occurs in less than 0.5 Gyr (see their figure 
7). This timescale is still younger than the typical estimated 
ages of high redshift compact ETGs. However, it should be 
pointed out that the details of this result depend also on the 
adopted recipes for stellar lifetimes and yields. These ingre- 
dients have some uncertainty (e.g. Romano et al. 2005). As 
a result, it cannot be totally excluded that passively evolv- 
ing ETGs lost 20-30% of the residual baryonic mass even 
~ 0.5 Gyr after the end of their main star forming phase. 
This would produce a small (due also to the relatively inef- 
ficient size increasing effect of slow mass loss, Figs. [1] and [2]), 
but not negligible, contribution to the claimed size evolu- 
tion. Moreover, in this case the uncertainty on the relatively 
difficult estimation of ages could have some importance, at 
least for the youngest observed high-z compact ETGs. 

To better illustrate the above points, we show in Fig. [7] 
the result of a sample simulation, applied to a specific semi- 
analytic galaxy formation model including both processes, 
namely QSO driven galactic wind and mass loss from stars 
due to stellar evolution. The figure displays the time evo- 
lution of baryonic mass as predicted by the G04 spheroid- 
SMBH co-evolution model (with the parameters as in Lapi 
et al. 2006), in a 1O 13 M DM halo that virializes at z = 4, 
together with the corresponding increase in size, computed 
with the procedure described in the present paper. The 
abrupt decrease of the mass after ~ 0.3 Gyr marks the 
ejection of gas by the AGN-driven wind. Note that, though 
this is a relatively fast process, it still does occur, according 



to the adopted recipes, on a timescales of a few tdya- This 
holds true over a wide range of model parameters. As a con- 
sequence, the corresponding puffmg-up is milder (a factor 
~ 1.5) than the maximal one (for a given e), achieved when 
At = 0. The subsequent slow decrease of mass and moderate 
increase in size (a further factor ~ 1.35, for a total expansion 
of ~ 2), is due to stellar mass returned to the ISM, under 
the assumption that the galaxy potential cannot retain it. 
The size expansion achieved after the epoch in which stellar 
populations are older than ~ 0.5 — 1 Gyr is 25% - 20%, 

In conclusion, the putative puffing up related to large 
scale galactic winds, quickly ejecting a substantial fraction 
of baryonic mass, can be an important phenomenon, but 
is still not observed. In particular, it cannot be invoked to 
explain the size evolution of ETG from z ~ 2.5 to z = 
observed in the presently available data. By converse, the 
secular adiabatic expansion, related to the mass returned 
to the ISM by stars during the final stages of their evolu- 
tion, could contribute, but not dominate, the observed size 
evolution of ETGs. Nevertheless, it is relevant to further in- 
vestigate also this contribution, since it seems that none of 
the processes or biases considered so far can explain alone 
this evolution (e.g. Hopkins et al. 2010). 



ACKNOWLEDGMENTS 

C.R-F. and G.L.G. acknowledge warm hospitality by INAF- 
Trieste and IATE-Cordoba, respectively, during the devel- 
opment of the present work. 

This work has been partially supported by the Consejo 
de Investigaciones Cientfficas y Tecnicas de la Republica 
Argentina (CONICET) and the Secretana de Ciencia y 
Tecnica de la Universidad Nacional de Cordoba (SeCyT). 

We thank Giuseppe Murante for carefully reading the 
manuscript and for useful suggestions, and the anonymous 
referee for proposing several improvements of the paper. 

REFERENCES 

Baumgardt, H., & Kroupa, P. 2007, MNRAS, 380, 1589 
Benson A. J., Bower R. G., Frenk C. S., Lacey C. G., Baugh 

C. M., Cole S., 2003, ApJ, 599, 38 
Benson A. J., 2010, PhR, 495, 33 
Bezanson, R., et al. 2009, ApJ, 697, 1290 
Biermann, P., & Shapiro, S.L. 1979, ApJ, 230, L33 
Boily, CM., & Kroupa, P. 2003, MNRAS, 338, 673 
Boylan-Kolchin, M., Ma, C.-R, & Quataert, E. 2006, MNRAS, 

369, 1081 

Buitrago, F., et al. 2008, ApJ, 687, L61 

Cassata P., et al., 2010, ApJ, 714, L79 

Cattaneo, A., et al. 2006, MNRAS, 370, 1651 

Cimatti, A., et al. 2008, A&Ap, 482, 21 

Ciotti L., Ostriker J. P., Proga D., 2009, ApJ, 699, 89 

Ciotti L., van Albada T. S., 2001, ApJ, 552, L13 

Daddi E., et al., 2005, ApJ, 626, 680 

Damjanov, I., et al. 2009, ApJ, 695, 101 

Fabian A. C, 1999, MNRAS, 308, L39 

Fan, L., Lapi, A., De Zotti, G., & Danese, L. 2008, ApJ, 689, 
L101 

Fan, L., Lapi, A., Bernardi, M., Bressan, A., Dc Zotti, G., & 

Danese, L. 2010, ApJ, 718, 1460 
Geyer, M.P., & Burkert, A. 2001, MNRAS, 323, 988 



10 Ragone-Figueroa & Granato 



Goodwin, S.P. & Bastian, N. 2006, MNRAS, 373, 752 

Granato G. L., Silva L., Monaco P., Panuzzo P., Salucci P., De 

Zotti G., Danese L., 2001, MNRAS, 324, 757 
Granato, G.L., et al. 2004, ApJ, 600, 580 (G04) 
Haring, N., & Rix, H.-W. 2004, ApJ, 604, L89 
Hernquist L., 1990, ApJ, 356, 359 
Hills, J.G. 1980, ApJ, 235, 986 

Hopkins P. F., Hernquist L., Cox T. J., Keres D., Wuyts S., 2009, 

ApJ, 691, 1424 
Hopkins, P.F., et al. 2010, MNRAS, 401, 1099 
Johansson, PH., Naab, T., & Burkert, A. 2009, ApJ, 690, 802 
Khochfar S., Silk J., 2006, ApJ, 648, L21 

Klypin A., Trujillo-Gomez S., Primack J., 2010, arXiv, 

larXiv:1002.3660l 
Komatsu, E., et al. 2009, ApJs, 180, 330 

La Barbcra F., dc Carvalho R. R., de la Rosa I. G., Sorrentino 

G., Gal R. R., Kohl-Moreira J. L., 2009, AJ, 137, 3942 
Lapi, A., et al. 2006, ApJ, 650, 42 
Lfpari, S., et al. 2009, MNRAS, 398, 658 
Longhetti, M., et al. 2007, MNRAS, 374, 614 
Mancini, C, et al. 2010, MNRAS, 401, 933 

McGrath E. J., Stockton A., Canalizo G., Iye M., Maihara T., 

2008, ApJ, 682, 303 
Monaco, P., Fontanot, F., & Taffoni, G. 2007, MNRAS, 375, 1189 
Naab, T., Johansson, P.H., & Ostriker, J.P. 2009, ApJ, 699, L178 
Navarro, J.F., Frenk, C.S., & White, S.D.M. 1997, ApJ, 490, 493 
Nipoti C, Londrillo P., Ciotti L., 2003, MNRAS, 342, 501 
Onodera M., et al., 2010, ApJ, 715, L6 

Oser L., Ostriker J. P., Naab T., Johansson P. H., Burkert A., 

2010, ApJ, 725, 2312 
Pipino A., Silk J., Matteucci F., 2009, MNRAS, 392, 475 
Prochaska J. X., Hennawi J. F., 2009, ApJ, 690, 1558 
Ryan R. E., Jr., et al., 2010, arXiv. larXiv: 1007. 14601 
Saracco P., Longhetti M., Andreon S., 2009, MNRAS, 392, 718 
Saracco P., Longhetti M., Gargiulo A., 2010, MNRAS, 408, L21 
Shen, S., et al. 2003, MNRAS, 343, 978 

Sijacki D., Springel V., Di Matteo T., Hernquist L., 2007, MN- 
RAS, 380, 877 
Silk, J., k, Rees, M.J. 1998, A&A, 331, LI 
Somerville, R.S., et al. 2008, MNRAS, 391, 481 
Springel V., 2005, MNRAS, 364, 1105 
Toft S., et al., 2007, ApJ, 671, 285 
Trujillo I., et al., 2006, ApJ, 650, 18 
Trujillo, I., et al. 2007, MNRAS, 382, 109 

Trujillo I., Ferreras I., de la Rosa I. G., 2011, arXiv, 

larXiv:1102.3398l 
Tutukov A. V., 1978, A&A, 70, 57 
Valentinuzzi T., et al., 2010, ApJ, 712, 226 
van der Wei, A., et al. 2008, ApJ, 688, 48 
van Dokkum, P.G., et al. 2008, ApJ, 677, L5 

Zhao, D.H., Mo, H.J., Jing, Y.P., & Borner, G. 2003, MNRAS, 
339, 12 

Zirm, A.W., et al. 2007, ApJ, 656, 66 



